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We derive the path integral of the semiclassical, one dimensional anharmonic Holstein model 
assuming that the electron motion takes place in a bath of non linear oscillators with quartic on 
site hard (and soft) potentials. The interplay between e-ph coupling and anharmonic force con- 
stant is analysed both in the adiabatic and antiadiabatic regime. In the latter we find much larger 
anharmonic features on the thermodynamic properties of low energy oscillators. Soft on site poten- 
tials generate attractive centres at large amplitude oscillator paths and contribute to the anomalous 
shape of the heat capacity over temperature ratio in the intermediate to low T range. This anhar- 
monic lattice effect is superimposed to the purely electronic contribution associated to a temperature 
dependent hopping with variable range inducing local disorder in the system. 
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I. INTRODUCTION 



There is at present a large interest on the effects of 
a strong electron-phonon coupling in a number of sys- 
tems ranging from dimer molecular junctions [1] to car- 
bon nanotubes [2] , from organic molecular crystals [3] , to 
DNA [4] and cuprate superconductors [5,6]. Several the- 
oretical studies have focussed on the interplay between 
e-ph coupling and non linearities in the framework of 
the Holstein model [7,8] investigating the phase diagram 
both in the adiabatic [9] and the antiadiabatic regime. 
Mainly in the latter anharmonic effects are believed to 
be large [10] thus offering a picture to explain the high 
Tc of binary alloys such as superconducting MgB2 with 
a small Fermi energy and sizeable e-ph coupling [11]. 

The path integral formalism provides a powerful 
method to study quantum systems in which a particle 
is non linearly coupled to the environment [12-14]. A 
previous path integral analysis [15] has pointed out how 
the phonon dispersion, that has to be taken into account 
in the computation of the ground state properties of the 
Holstein Hamiltonian [16-18], induces non local e-ph cor- 
relations which renormalize downwards the effective cou- 
pling and ultimately broaden the size of the polaronic 
quasiparticle. This explains why the polaron mass in a 
dispersive Holstein model [19] is lighter than in a dis- 
persionless model [20] . Also the thermodynamics of the 
Holstein Hamiltonian can be computed within a disper- 
sive model which accounts for the lattice structure [15]. 

The Holstein diatomic molecular model was originally 
cast [21] in the form of a discrete non linear Schrodinger 
equation for electrons whose probability amplitude at a 
molecular site depends on the interatomic vibration co- 
ordinates. The non linearities are tuned by the e-ph cou- 
pling [22] whose strength drives the crossover between a 
large and a small polaron state according to the degree 
of adiabaticity and the dimensionality of the system [23] . 
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In the Holstein model, the phonon thermodynamics is 
not affected by e-ph induced anharmonicities [24]. This 
follows from the fact that the Holstein perturbing source 
current is local in time and it does not depend on the 
electron path coordinate. As a consequence, in the total 
partition function, electron and lattice degrees of freedom 
are disentangled and the latter can be integrated out an- 
alytically as far as a harmonic lattice model is assumed. 
However non linearities may arise in Holstein like systems 
also by virtue of on site potentials dependent on the lat- 
tice structure and, in principle, independent of the e-ph 
coupling. We are thus led to investigate the thermody- 
namics of the anharmonic Holstein model with a quartic 
on site potential which may be repulsive or partly attrac- 
tive according to the sign of the force constant and the 
amplitude of the lattice displacement paths. The path 
integral approach permits to monitor the physical prop- 
erties for any value of the coupling strengths [25]. The 
presence of a 0^ potential may in turn affect also the e-ph 
interactions and sinergically interfere on the equilibrium 
properties of the system. This is the focus of the present 
paper. Section II presents the Hamiltonian model while 
the path integral method is briefly described in Section 
III. The derivation of the total partition function of the 
system is presented in Section IV and Section V contains 
the discussion of the physical results. The conclusions 
are drawn in Section VI. 



II. THE ANHARMONIC HOLSTEIN MODEL 

We consider the one dimensional anharmonic Holstein 
Hamiltonian consisting of: i) one electron hopping term, 
ii) an interaction which couples the electronic density 
(///i) to the lattice displacement ui at the 1-site, iii) a 
bath of N identical dispersionless anharmonic oscillators 
with mass M and frequency co : 

<l,m> 
N 

1=1 

1=1 1=1 

g = gV2Muj (0.1) 

the sum <, > is over z nearest neighbors and t is the 
tight binding overlap integral, g is the e-ph coupling 
in units of Huj. Choosing the atomic mass M of or- 
der 10^ times the electron mass, we get g ~ 1.1456 x 
gV^hcj [meVA~^] where Tiu is given in meV. 6, in 
units 'm,eVA~'^, controls the strength of the non lin- 
earities and determines whether the on site potential 
V{ui) = Moj'^uf/2 + Suf/A is hard {5 > 0) or soft 
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{S < 0). In the latter case, V{ui) attains the maxi- 
mum at = M(jj'^/\d\ and the inflection point occurs 
at uf = Mlo'^/3\5\. The condition \5\uf > 2Muj^ yields 
an attractive on site potential. Then, the range of the 
atomic path amplitudes generating attractive scattering 
centres depends on the value of the anharmonic force con- 
stant. For \5\ > 2Mu!^ the potential becomes attractive 
for a portion of large amplitude atomic paths while small 
amplitude paths weigh the repulsive range. 

In the following computation of the electron path in- 
tegral coupled to the anharmonic oscillator, after setting 
the potential parameters, we select at any temperature 
the class of atomic paths which mainly contribute to the 
euclidean action. As the distribution of the path ampli- 
tudes has a cutoff on the scale of the lattice constant, say 
i.e. uf < lA^, the on site potentials are always bound 
from below also in the attractive cases. 



III. THE PATH INTEGRAL METHOD 

The Holstein Hamiltonian in (1) can be mapped onto 

the time scale according to space-time mapping tech- 
niques extensively described in previous works [15,26,27] 
and hereafter outlined. 

Defining x(t) and ^(t') as the electron coordinates at 
the I and m lattice sites respectively, if in (1) transforms 
into 

H^{t,t') = -t{p{x{T))f{y{r')) + p{y{T'))f{x{T))) 

(0.2) 

where r and r' are continuous variables (e [0, /3]) , with 
P being the inverse temperature. After setting r' = 0, 
y(0) = and taking the thermal averages for the electron 
operators over the ground state of the Hamiltonian one 
gets the average electron hopping energy per lattice site: 

h%x{T)) ^ ^^^5^ = -t{G[-xiT),-T]+G[xir),T]) 

(0.3) 

where, G[x{t),t] is the electron propagator at finite 
temperature. 

By treating the lattice displacements in (1) as r— de- 
pendent classical variables, ui — » u{t), we obtain from 
jje-ph ^j^^ ^jjg averaged e-ph energy per lattice site 
which is identified as the perturbing source current j(r) 
in the path integral method: 

i(T) ^ ^ ^^^^ = gu{T) (0.4) 

As the Hamiltonian model assumes a set of identical 
oscillators we study the path integral for the electron 
coupled to a single anharmonic oscillator of the bath. 
The path integral reads: 
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< x{/3)\x{0) >= j Dx{T)exp - dTE{x{T)) 
j Du{T)exp — J dTO{u{T)) 



E{x{T))='-fx\T)+h%x{T)) 

0{u{t)) = ^(u\r)+M{T)) + ^-u\t) + j(r) 

(0.5) 

where the kinetic term (me is the electron mass) is 
normalized by the functional measure of integration over 
the electron paths. 

Since the electron hopping does not induce a shift of 
the oscillator coordinate the Holstcin e-ph interactions 
are local in time and, in the semiclassical treatment, 
the source current j(t) is independent of the electron 
path. As a consequence oscillator and electron coordi- 
nates appear disentangled in (5) while the coupling oc- 
curs through the parameter g. 



IV. THE PARTITION FUNCTION 



The quantum statistical partition function Zt is de- 
rived by integrating (5) after imposing periodicity condi- 
tions, /3 is the period, both on the electron and oscillator 
paths: 

Zt= j dx< x{f5)\x{0) >= Za X Zosc 



Zei = j> Dx{T)exp 



dTE{x{T)) 



^osc — 



(0.6) 



where § Dx{t) and § Du{t) are the functional mea- 
sures of integration. 

The electronic contribution is computed by ex- 
panding the paths in Fourier components 

Mp 

x{t) = Xo+^ [rm COs(yrnT) + Sm Sm{VmT)) 

m— 1 

= 2TO7r//3 



(0.7) 



and taking the following measure of integration 

a/2 



j Dx{t) 



(2A„J(2A^-+i) 



/ 

J —t 



dxo 
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/.OO /'OD 

{2-KmY i drm / ds 



(0.8) 



The cutoffs over the Fourier coefficient integrations 
have to ensure proper normaUzation of the kinetic term 
in absence of hopping processes. Thus Zf,i transforms 
into 



- j(.M..r) \_^^^ n ds. 



_ ^/2 r^/^ 
(2Ato 

3 Mf 



X exp --2- ^ m2(r^ + Sm) - / dr/i^(a:(r)) 

(0.9) 

with A oc A„j^ [27] indicating that large amplitude elec- 
tron paths have to be selected at low temperatures where 
the quantum effects are larger. Two Fourier components, 
Mp = 2, suffice to attain stable results as h'^{x{T)) de- 
pends smoothly on the electron path. The hopping term 
accounts for the deviation from the Gaussian behavior. 
Numerical analysis shows that, for any choice of path pa- 
rameters, h'^{x{T)) decreases by decreasing the tempera- 
ture but its overall contribution to the electron action is 
substantial also at low T. 

Let's focus now on the anharmonic oscillator term Zq^q. 
The oscillator path is expanded in Np Fourier compo- 
nents 



Nf 

E 

n=l 



i{t) = Mo + ^ {an cos{u)„t) + b„ sin(a;„T)) 

(0.10) 



with Matsubara frequencies Un ~ 2mT / (3 and coeffi- 
cients a„ = 3?u„, bn = — satisfying the conditions 
a„ = a_„ and 6„ = —h-n- The latter are consistent 
with the choice of real paths and simplify the following 
T integration of the on site potential. 

Note that the periodicity property, u{t) = u{t + (3), 
would be fulfilled also taking the very a„ coefficients in 
(10) [28]. However such a choice would not permit to fit 
with accuracy the harmonic oscillator partition function 
which is known exactly: Zh = [2 sinh(/3a;/2)] . In- 
fact, in the path integral method, the harmonic partition 



function Zf^^ reads 



Z^^ = j> Du{T)exp 



(0.11) 

and taking the functional measure 
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Yj / dUo 

J —OO 

J'' poo poo 

JJ(27rn)^ / dan / dbn 

1 J — oo J —oo 



Nf 



n=l 



(0.12) 



with Am = ^J^^hFp/M, one gets from (10)-(12): 



7PI 



(2n7r)2 



(3uj n (2n7r)2 + (/3a;)2 



(0.13) 



Instead, dropping the 6„ terms in (10) and (12), one 
would get the square root of the product series in (13) 
which does not yield a reliable fit of even for large Np. 
Note that at high T, the condition 2n7r ^ /3cj is fulfilled 
for small integers n, hence the main contribution to Z^^ 
is given by the 1 / /Su factor which stems from the / duo in 
(12). This is consistent with the expectation that high T 
paths are well approximated by their /3-avcragcd value Uq 
whereas fluctuation effects become increasingly relevant 
towards the low T regime in which Np rapidly grows. 
Np clearly varies also with the oscillator energy while the 
shape of Np{lu,T) may differ according to the harmonic 
function {Z^, harmonic free energy or specific heat) one 
chooses to fit. 

The anharmonic partition function Zosc in (6) can be 
worked out analytically using (10), (12) and performing 
the time integration of the oscillator functional 0(m,(t)). 
This permits to get an insight into the role of the non 
linear terms. The lenghty calculation yields: 



V2 



^"'^ ~ (2Am)(2^' 

Np 



/oo 
dUo exp{-PgUo - nul - PSu^/A) 



n(27m)^ 



ra=l 



da„ exp 



Nf 



n=l 



Nf 



Nf 



l35Uo / \ , ^ ^( ^ 

-— — 2^ c(n,m) + — 2^ d{n,m,p) 



m,p=l 



dbn 



Nf 

E 

n=l 



[(7n + ^PSul/A)bl 



Nf 



Nf 



3f3Suo ^ g(^^^) _^ ^ {6f{n,m,p) - g{r 



hP)) 



m,p=l 



K = 7r(/3w) V2A^ 

7„= 7r((27rn)2 + (/3a;)2)/4Ai^ 

c(n, m) = anam{an+m + an-m. 

d{n,Tn,p) = anamap{an+m+p ^ 
e{n,m) = a„6„i(6„+m - 6„-m) 
f{n,m,p) 



^n—m-\-p 



H" ^p—n-\-m) 



hn- 



■m+p 



bp—n—m ~i~ bp—n-\-m'} 



g{n, m,p) = bnbmbp{bn-\-m-\-p - bn-m+p + bp- 



^p—n-\-m J 



6 



(0.14) 



In c{n, m) the coefficients aj {j = n + m,, n — m) vanish 
if: j < or j > Np. In d{n, m,p), Uk ^ {k = n + m + 
p, n — m + p, p — n — m, p — n + m) '^=> 1 < A; < Np- 
Analogous conditions hold for the coefficients bj (/;) in the 
e{n,m), f(n,m,p) and g{n,m,p) functions. 

Note that the effective e-ph coupling g is associated 
only to the r-independent component Ug, that is to the 
/3-averaged displacement path [15]. This follows from the 
fact that the perturbing source current is non retarded in 
the Holstein model as a consequence of the local nature 
of the e-ph interaction. 

The quartic potential induces a strong mixing of the 
Fourier components of the path that highly complicates 
the numerical problem. Thus the value of Np appears to 
be crucial in the computation. We determine Np{lu,T) 
by fitting (with an accuracy of 2 • 10~^) the exact har- 
monic free energy {Fh = — ln{Zh)/f3), through the path 
integral harmonic free energy _F)f^ obtained from (13). 
As an example, for the oscillator with u = 20meV, 
Nf{T = lOK) = 59 and Np{T ^ 200K) = 8. 

Inspection of (14) offers the key to perform reliable 
path integrations according to the sign of the 0^ poten- 
tial. At high temperature, a large contribution to Z^g^ is 
expected to come from the paths having Ug which maxi- 
mizes exp(K/(Mo)) with f(uo) = —{aUo + Ug + bUg), a = 
Pg / n\ b = [3S/Ak. In general, we find that for a 
hard (soft) potential, the dwo-integration has to be car- 
ried out along the < {uq > 0) axis, with cutoff 
|uo|ma:r ^ 0.6/\/k. This permits to include the set of 
paths which mainly contribute to the euclidean oscillator 
action. On the Fourier coefficients integrals, / rfa„ J c?6„, 
we set the cutoffs \an\max, \bn\max ^ O.Q/^/j^ both for 
hard and soft potentials thus achieving numerical conver- 
gence and correct computation of the Gaussian integrals 
once the non linearities are switched off. It turns out, see 
the definitions in (14), that the cutoffs on the oscillator 
path integration are increasing functions of temperature 
(oc a/T) consistently with the physical expectations of 
large amplitude displacements at high T. As the path 
displacements encounter an upper limit due to the cut- 
offs, u{t) < \uo\max + "^Y^n^i Wn\max, the distribution 
of on site potentials has a lower limit also in the case of 
soft and attractive non linearities. This avoids numerical 
divergences and makes the problem physically meaning- 
ful. 



V. RESULTS 

We test the relevance of the non linearities on the equi- 
librium thermodynamics of the system and present the 
calculation for the heat capacity in the intermediate to 
low temperature range. 

Figures 1 show the behavior of a low energy {u = 
20meV) oscillator without {g = 0) and with {g = 
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2, 4) coupling to the electronic subsystem in the adia- 
batic regime: t/tj = 5. Figs. 1(a), (b) assume an an- 
harmonic potential with positive quartic force constant 
6 = 10^, lO'^meyA"'* respectively. Also the harmonic 
heat capacity is reported on for comparison. The hard 
potential lifts the free energy over the harmonic values 
with a more pronounced enhancement at increasing T 
and for larger 6. The effect on the free energy second 
derivative is however scarce and essentially consists, see 
Fig. 1(a), in a slight increase ( decrease) of the heat ca- 
pacity at low (high) temperatures. The reduction of the 
heat capacity (with respect to the harmonic plot) at in- 
termediate and high T is more evident in Fig. 1(b) where 
the quartic force constant is larger. This result is in ac- 
cordance with diagrammatic perturbative treatments of 
the anharmonic crystals which predict negative contri- 
butions to the constant volume specific heat arising from 
positive force constants in the quartic potential [29,30]. 
The e-ph coupling strength also tends to decrease the os- 
cillator heat capacity but the overall effect is small as 
the insets in the figures show: infact the electronic term 
dominates the total heat capacity C and the oscillator 
contribution is not distinguishable in the plots of the C 
over temperature ratios versus T. The low temperature 
upturn is due to the large electron energy term in (9) and 
precisely ascribable to the feature of the variable range 
(on the r scale) of the electron hopping, captured by the 
path integral formalism [27]. 

The cases of a soft on site potential are reported on in 
Fig. 1(c) with S = -lO^meVA-^ and Fig. 1(d) with 5 = 
— 10'^meVA~^. The characteristic potential parameter is 
Mw^/I^l ~ 0.5 and O.OSA^, respectively. The potential 
V{u{t)) = Ma;2M(T)2/2-h(5M(T)^/4 is attractive for those 
paths such that u{t)'^ > 2Mtj'^/\S\. 

At any T, we integrate over a distribution of time de- 
pendent potentials. Thus, at a given lattice site, the 
electron may experience an attractive or repulsive scat- 
tering centre according to the size of the atomic path. 
As an example, at T = 200^, we get a maximum path 
Umax such that u^^j, ^ 0.15A^. This guarantees that 
V{u{t)) is generally repulsive in Fig. 1(c) and attractive 
for a broad class of paths in Fig. 1(d). 

The effects on the oscillator heat capacity are twofold: 
i) the soft potential enhances the heat capacity mainly in 
the low T range with respect to the hard potential and 
this feature is much more pronounced in Fig. 1(d); ii) the 
trend of the e-ph coupling is opposite to that observed 
in Figs. 1(a), (b): now by increasing the g values one gets 
higher heat capacities although the size of this effect is 
small on the scale of the electronic terms in the adiabatic 
regime as revealed by the C/T plots in the insets. 

Let's come to the antiadiabatic regime {t/uj = 0.5) 
discussed in Figures 2, where a low harmonic energy 
{lo = IQmeV) is assumed to emphasize the size of the an- 
harmonicity together with a very narrow electron band. 
A hard on site potential with 5 = lO^raeV A~'^ is taken in 
Fig. 2(a): the shape of the oscillator heat capacity signals 
the effects of the non linearities which flatten the curve at 
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intermediate T and enhance the heat capacity also at low 
T with respect to the corresponding case of Fig. 1(a). As 
the Debye temperature is now smaller (than in Fig. 1(a)) 
the hard anharmonicity decreases the constant volume 
heat capacity with respect to the harmonic plot over a 
broader temperature range. 

The e-ph coupling plays a minor role also in antiadia- 
batic conditions. The anharmonic contribution is visible 
in the total heat capacity as the inset makes evident al- 
though the dominant electronic feature persists in the 
low T upturn of C/T. In Fig. 2(b), we assume g = 2 
and consider two cases of soft potential: the oscillator 
anharmonicity becomes relevant and such to modify the 
shape of the anomalous upturn in the total heat capac- 
ity. An enhancement of the C/T values is observed at 
intermediate and low T and, in the case of the largest \6\ 
generating a soft attractive potential, the oscillator heat 
capacity Cqsc yields an upturn in Cqsc/T independently 
of the electronic term. 



VI. CONCLUSIONS 

We have studied the path integral of the one dimen- 
sional non linear Holstein model in which a set of disper- 
sionlcss oscillators provides the environment for the elec- 
tron. The model is semiclassical as the lattice displace- 
ments are treated classically while the electron operators 
are thermally averaged over the ground state Hamilto- 
nian. The e-ph coupling of the model is local and gener- 
ates a perturbing current which linearly depends on the 
oscillator path amplitude u(r) where r is the time (or 
inverse temperature) of the Matsubara Green functions 
formalism. The anharmonicity on the lattice site is mod- 
elled through a (j)^ potential that may result attractive 
for a set of displacement paths in the case of a negative 
quartic force constant (soft potential). We have derived 
the path integral of the interacting system and computed 
the total partition function selecting, as a function of the 
temperature (T < 200-/^), both the electron and oscil- 
lator paths which yield the largest contribution to the 
action. While quantum electron paths have increasing 
amplitudes at decreasing temperatures, the atomic dis- 
placements are growing functions of T. This relevant 
physical feature is accounted for in our model as the cut- 
offs on the electron path integration are proportional to 
the electron thermal wavelength whereas, on the atomic 
path integration, we find cutoffs proportional to ^/T. 

The oscillator partition function includes the effect of 
the coupling to the electron subsystem but the on site 
anharmonicitics play a major role mainly when the po- 
tential is soft and the harmonic energy is low. Among the 
thermodynamic properties we have chosen to present the 
heat capacity C in view of the upturn in the C/T behav- 
ior due to the low T electron hopping tuned by the value 
of the overlap integral in the Hamiltonian model. The 
computation is highly time consuming especially because 
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of the strong mixing of the path Fourier components gen- 
erated by the non hncar potential. 

In general we find that: i) in the case of a hard quartic 
potential, switching on the e-ph coupling leads to lower 
the free energy and its second derivative, ii) when the 
quartic potential is soft, e-p/i coupling and anharmonicity 
act sinergically enhancing the thermodynamic functions. 

These results can be understood on general physical 
grounds. Infact, a hard quartic potential shifts the char- 
acteristic phonon frequency upwards thus hardening the 
spectrum and broadening the size of the quasiparticle. If 
the e-ph coupling is enhanced (independently of the on 
site anharmonicity) the magnitude of the source which 
causes the lattice distortion becomes larger. This further 
hardens the vibrations and leads to decrease the anhar- 
monic heat capacity over the whole temperature range. 
Instead, when the quartic potential is soft the phonon 
frequency is lowered and the oscillator potential well is 
more flexible. In this case, larger e-ph coupling strengths 
favour the self-trapping of quasiparticles with heavier ef- 
fective masses. This is physically equivalent to soften the 
phonon spectrum and enhance the heat capacity. 

The electron contribution to the heat capacity is dom- 
inant in the adiabatic regime whereas antiadiabatic sys- 
tems are expected to present significant anharmonic cor- 
rections. Infact, in the antiadiabatic regime the quasi- 
particle is a small size object on the lattice scale and 
the electron energy associated with the overlap integral 
is small. Thus, this regime proposes a physical picture in 
which the electron hardly hops from site to site and its 
effective mass becomes heavier on the scale of the atomic 
mass. But a potential well generated by "lighter oscilla- 
tors" is more sensitive to on site anharmonic effects. 

In particular, soft potentials increase the heat capaci- 
ties over the harmonic values and reinforce the upturn in 
the C/T versus T plots when the on site anharmonicity 
is such to produce attractive potentials for a set of lat- 
tice displacement paths. Since the path amplitudes are 
larger at high T, soft attractive potentials induce rapidly 
increasing phonon heat capacities at growing T as shown 
in Fig. 2(b). 

Thus our path integral investigation and the thermo- 
dynamical results point to a complex role of the lattice 
anharmonicities in the one dimensional Holstein model 
and suggest that on site potentials may be experienced 
as attractive or repulsive according to the temperature 
and the amplitude of the atomic path. Such potentials 
may provide scattering centres generating local disorder 
whose effect on the system thermodynamics is superim- 
posed to the disorder induced by the hopping of electrons 
with variable range. 
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FIG. 1. (Color online) Anharmonic Oscillator Heat Capac- 
ities versus temperature for three values of e-ph coupling g and 

oscillator energy a; — 20meV. (a) Hard potential with force 
constant 5 = lO^meVA-'^; (b) 5 = lO^meVA"''; (c) soft po- 
tential with 6 = -lO^meVA-'^; (d) 5 = -W^meVA-*. The 
harmonic heat capacity is plotted in (a) and (b) for compar- 
ison. The insets show the Total (electronic plus anharmonic 
oscillator) Heat Capacity over temperature ratios in the adia- 
batic regime t/uj = 5. 
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FIG. 2. (Color online) Anharmonic Oscillator Heat Capac- 
ities versus temperature in the antiadiabatic regime tju) = 0.5 
and oscillator energy u} = WmeV. The harmonic heat ca- 
pacity is also plotted, (a) Hard potential force constant 
6 = l0^meVA~'^ with three values of e-ph coupling g. (b) 
Two soft potential force constants at fixed e-ph coupling. The 
insets show the Total (electronic plus anharmonic oscillator) 
Heat Capacity over temperature ratios. The electronic contri- 
bution is plotted separately for comparison. 
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